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ABSTRACT 

two  dimensional  elastodynamic  problem  of  a  semi-infinite  plate  containing  an  edge 
crack  is  considered.  Initially,  the  plate  is  stress  free  and  at  rest.  To  simulate  the  asymmetric 
impact  of  a  projectile  on  the  cracked  edge  of  the  plate,  a  normal  velocity  is  suddenly  imposed 
on  the  boundary  of  the  plate  on  one  side  of  the  edge  crack.  The  boundary  of  the  plate  and  the 
crack  faces  are  otherwise  traction  free.  Due  to  the  nature  of  the  loading,  a  combined  mode  I 
and  mode  II  transient  stress  intensity  factor  is  induced  at  the  crack  tip.  This  stress  intensity 
factor  history  is  determined  exactly  by  linear  superposition  of  several  more  readily  obtainable 
stress  wave  propagation  solutions,  including  a  fundamental  solution  arising  from  a  particular 
problem  in  the  dynamic  theory  of  the  elastic  dislocations.  The  stress  intensity  factor  histories 
are  determined  for  the  time  interval  from  initial  loading  until  the  first  scattered  wave  at  the 
crack  tip  is  reflected  at  the  edge  and  arrives  at  the  crack  tip  again.  For  an  applied  velocity  with 
step  function  time  dependence,  the  mode  I  stress  intensity  factor  monotonically  decreases  and 
that  of  mode  II  increases  until  the  arrival  of  the  Rayleigh  wave,  it  then  slightly  decreases  for 
short  duration,  and  then  increases  again.  In  experiments  on  fracture  initiation  in  a  high  strength 
steeljbased  on  essentially  this  specimen  and  loading  configuration,  KalthofT  and  Winkler  (1987) 
reported  tharthe  fracture  cither  grew  off  at  an  angle  as  a  tensile  crack  or  strain  ahead  as  a  shear 
fracture,^depending  on  the  intensity  of  the  applied  loading.  The  observations  arc  considered  in 
light  of  the  solution  reported  hcre^'Yt*!;?/ 


1.  INTRODUCTION 


An  experimental  technique  has  been  proposed  by  Kalthoff  and  Winkler  (1987)  and  Kalthoff 
(1987)  for  subjecting  edge  cracks  in  plate  specimens  to  very  high  rates  of  loading  that  result 
in  a  crack  tip  deformation  field  that  is  predominantly  mode  II,  that  is,  the  in-plane  shearing 
mode.  A  plate  specimen  with  two  parallel  edge  cracks  or  notches  is  impacted  by  a  cylindrical 
projectile  of  diameter  equal  to  the  spacing  between  the  cracks,  as  shown  schematically  in  Fig.  1. 
Upon  impact,  the  projectile  produces  a  compressive  wave  in  the  part  of  the  specimen  between 
the  two  cracks,  and  this  wave  propagates  toward  the  crack  tips  while  merely  grazing  the  two 
crack  faces  that  border  this  part  of  the  specimen.  Upon  arrival  of  the  wave  at  the  crack  tips, 
the  stiffness  constraint  of  the  remainder  of  the  specimen  is  encountered  and  relatively  large 
shear  stress  is  induced  on  the  crack  planes  ahead  of  the  crack  tips  resulting  in  transient  mode 
n  stress  intensity  factors.  The  same  idea  has  been  applied  by  Kalthoff  and  co-workers  for  the 
case  of  a  plate  with  a  single  edge  crack. 

The  wave  motion  in  the  part  of  the  specimen  between  the  cracks  upon  projectile  impact  is 
actually  much  more  complicated  than  the  above  discussion  would  suggest.  A  plane  compressive 
wave  is  indeed  generated  upon  impact.  Because  the  crack  faces  are  free  of  traction,  cylindrical 
unloading  waves  immediately  begin  to  propagate  outward  from  tlie  corners  where  the  crack 
faces  intersect  the  impacted  edge  of  the  plate.  The  wavefront  diagram  is  shown  in  Fig.  2.  These 
unloading  waves  have  two  main  effects.  They  tend  to  erode  the  strength  of  the  compressive 
plane  wave  as  it  propagates  along  the  crack  faces  and  they  tend  to  result  in  outward  bulging  of 
the  crack  faces  as  the  intense  compression  of  the  plane  wave  is  relieved  due  to  the  presence  of 
the  free  surface.  These  features  are  evident  in  some  early  work  by  Skalak  (1957)  on  longitudinal 
wave  propagation  in  a  plate  due  to  impact  of  the  edge  of  the  plate.  The  net  result  of  the  impact 
loading  is  to  induce  a  transient,  mixed-mode  stress  intensity  factor  field  at  the  crack  tip.  The 
mode  II  stress  intensity  has  a  large  component  due  to  the  main  compressive  wave,  but  a  time 
dependent  mode  I  stress  intensity  factor  is  also  generated.  In  view  of  the  tendency  for  the  crack 
to  close  noted  above,  the  mode  I  stress  intensity  factor  will  likely  be  negative. 

Data  are  reported  by  Kalthoff  (1987)  for  experiments  on  specimens  of  a  high  strength 
maraging  steel.  The  plate  dimensions  were  approximately  1  cm  X  10  cm  X  20  cm,  and  the  crack 
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length,  crack  separation  distance  and  projectile  diameter  were  all  about  5  cm.  Impact  velocities 
were  in  the  range  from  about  10  m/s  to  100  m/s.  A  particularly  interesting  feature  of  the 
data  is  that  two  completely  different  failure  modes  were  observed  to  occur,  depending  on  the 
impacting  speed  of  the  projectile.  For  the  lowest  impact  velocities,  fracture  occurred  in  the 
form  of  planar  crack  growth  in  a  direction  inclined  at  about  70°  to  the  original  crack  plane.  For 
impact  speed  greater  than  about  20  m/s,  on  the  other  hand,  failure  occurred  by  the  crack-like 
growth  of  a  shear  band  in  a  direction  almost  coincident  with  the  original  crack  plane. 

The  purpose  here  is  to  compute  the  transient  mixed  mode  stress  intensity  factors  through 
analysis  of  a  boundary  value  problem  intended  to  model  the  experimental  conditions  in  order 
to  provide  a  broader  basis  on  which  to  interpret  the  experimental  results.  The  analysis  here  is 
based  on  the  assumption  of  elastic  material  response,  and  study  of  the  influence  of  crack  tip 
plasticity  is  underway. 

The  particular  edge  cracked  plate  problem  to  be  analyzed  is  shown  schematically  in  Fig.  3. 
A  rectangular  i,  y-coordinate  system  is  introduced  in  the  plane  of  the  plate,  and  the  plate 
occupies  the  half  plane  x  >  -1.  A  crack  with  traction  free  faces  extends  normally  inward  a 
distance  I  from  the  edge  of  the  plate.  Thus,  the  crack  tip  is  at  the  origin  of  coordinates  as 
shown  in  Fig.  3.  The  normal  stress  component  (Txx  vanishes  on  the  edge  of  the  plate  for  y  >  0 
and  a  time  dependent  normal  velocity  is  imposed  on  the  edge  of  the  plate  for  y  <  0.  The  shear 
traction  is  assumed  to  vanish  everywhere  on  the  plate  edge. 

To  solve  the  problem,  the  longitudinal  displacement  potential  </>  and  the  shear  displacement 
potential  V’  are  introduced  through  the  Helmholtz  decomposition  of  the  displacement  vector. 
The  components  of  displacement  in  the  coordinate  directions  are  derived  from  the  potentials 
by  differentiation  according  to 


«!  =  </>,!  + t/’.v  My  =  0,y  -  (1-1) 

where  the  familiar  comma-subscript  notation  is  used  to  denote  partial  differentiation. 

Each  of  these  potential  functions  satisfies  a  two  dimensional  wave  equation, 

*P,xx  ^,yy  ~  ®  4^, it  —  0  V’.xar  "h  ti^.yy  ~  ^  4’, it  ~  0  (1'2) 
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where  a  =  l/c^  =  \/pl(X  +  2/i)  and  b  =  1/c,  =  y/pjp  are  the  inverse  dilatational  and  shear 
waves  speeds,  respectively,  in  terms  of  the  Lame  elastic  constants  A,  p.  and  the  mass  density  p. 
For  future  reference,  c  =  Ijcji  is  the  inverse  Rayleigh  wave  speed. 

The  stress  components  can  be  expressed  in  terms  of  the  displacement  potentials  by  means 
of  Hooke’s  law  for  the  material.  As  already  discussed  above,  the  portion  of  the  plate  edge  above 
the  crack  is  completely  traction  free.  On  the  other  hand,  the  portion  of  the  edge  below  the 
crack  is  free  of  shear  traction,  but  the  normal  speed  is  imposed  beginning  at  time  t  =  0.  Thus, 
the  solution  must  satisfy  the  boundary  conditions 

=  y>0 

=  0*  y<0  (1-3) 

<7xy(l,0*,t)  =  0,  — £  <  I  <  0 

where  v(t)  is  prescribed  for  t  >  0.  The  formulation  is  completed  by  specifying  zero  initial  data. 
Evidently,  the  stress  intensity  factors  will  be  identically  zero  until  time  t  =  tjcd- 

The  boundary  value  problem  formulated  above  is  linear.  Nonetheless,  it  is  sufficiently 
complex  so  that  known  direct  methods  of  solution  are  apparently  not  applicable.  By  studying 
the  features  of  a  sequence  of  sub-problems,  however,  it  will  be  shown  that  mathematically 
exact  expressions  for  the  mode  I  and  mode  II  stress  intensity  factors  can  be  found.  The 
property  of  linearity  is  exploited  by  superposition  of  these  solutions  in  a  certain  way.  The 
Laplace  transform  is  useful  in  approaching  some  of  these  problems.  The  dual  transform  of  the 
dilatational  displacement  potential  <f>{x,y,t)  o\z\  time  t  and  spatial  coordinate  x  is  defined  by 

/OO  ^OO 

(1.4) 

-OO  Jo 

where  the  transform  parameters  are  defined  to  be  s  and  s^,  respectively.  The  parameter  s  may 
be  viewed  as  a  sufficiently  large  positive  real  number  in  considering  transform  on  x. 

The  superposition  scheme  is  outlined  in  Figs.  4,  5  and  6.  In  Fig.  4,  it  is  shown  that  the 
asymmetric  impact  problem  can  be  viewed  as  a  superposition  of  a  symmetric  impact  problem 
(A)  and  an  antisymmetric  impact  problem  (B).  If  v{t)  is  tlic  imposed  normal  velocity  in  the 
problem  of  interest,  then  a  normal  edge  velocity  of  magnitude  |r(0  is  imposed  symmetrically 


«x(-£,y,£)=  /  v(T)dr, 
Jo 

<Txxix,0^,t)  =  0, 
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with  respect  to  the  crack  plane  in  (A)  and  antisymmetrically  with  respect  to  the  crack  plane 
in  (B).  Thus,  solutions  of  these  two  problems  is  required. 

Consider  the  symmetric  problem  (A).  It  is  shown  in  Fig.  5  that  this  problem  can  be 
decomposed  into  a  quarter  plane  problem  (C)  and  a  problem  concerned  with  the  motion  of 
dislocations  emanating  from  the  tip  of  a  crack  (D).  The  basic  idea  is  as  follows.  If  the  imposed 
velocity  begins  to  act  on  the  edge  x  =  —tof  the  quarter  plane  in  (C),  then  a  normal  displacement 
is  induced  on  the  edge  y  =  0.  However,  the  displacement  in  the  y-direction  on  y  =  0  for  the 
problem  of  interest  (A)  is  zero.  Consequently,  the  normal  displacement  induced  on  i  >  0,  y  =  0 
for  the  quarter  plane  problem  (C)  is  canceled  by  a  distribution  of  climbing  edge  dislocations 
emanating  from  the  crack  tip  in  (D)  with  a  net  displacement  distribution  on  y  =  0,  x  >  0  that 
exactly  negates  the  normal  displacement  in  (C).  The  use  of  dynamic  dislocation  solutions  in 
this  way  was  introduced  by  Freund  (1974)  in  order  to  obtain  stress  intensity  factor  solutions 
for  certain  problems  that  are  not  amenable  to  direct  analysis.  These  solutions  have  been  used 
to  interpret  experimental  data  by  Ravi-Chandar  and  Knauss  (1984)  and  Kim  (1985).  A  similar 
superposition  scheme  can  be  used  to  solve  the  antisymmetric  loading  problem. 

Finally,  the  quarter  plane  problem  can  be  solved  by  the  superposition  shown  schematically 
in  Fig.  6.  Case  (E)  is  the  problem  of  a  half  plane  occupying  x'  =  x  +  f  >  0  with  the  surface 
of  the  half  plane  subject  to  a  uniform  norma]  velocity  ^v{t).  A  result  of  this  loading  is  that  a 
normal  stress  is  induced  on  y  =  0  for  x'  <  cjt.  This  normal  stress  is  then  negated  by  considering 
the  problem  (F)  of  a  half  plane  occupying  y  >  0  subjected  to  a  surface  normal  traction  equal 
but  opposite  to  that  induced  on  y  =  0  in  (E).  This  approach  was  introduced  by  Wright  (1969) 
and  Freund  and  Phillips  (1969)  in  studies  of  grazing  incidence  of  a  stress  pulse  on  a  free  surface. 

2.  QUARTER  PLANE  PROBLEM 

In  this  section,  the  solution  of  the  quarter  plane  problem  designated  as  (C)  in  the  preceding 
discussion  is  determined.  The  plane  strain  deformation  of  an  elastic  material  occupying  the 
region  x'  >  0,  y  >  in  the  x,y-plane  is  considered,  as  indicated  in  Fig.  6.  The  prime  on  x  is 
dropped  for  the  development  in  this  section.  Initially,  the  material  is  stress  free  and  at  rest. 
At  time  /  =  0,  a  spatially  uniform  normal  displacement  is  imposed  on  the  edge  x  =  0,  y  >  0 
of  the  quarter  plane.  The  shear  traction  is  zero  on  this  edge,  and  the  traction  is  also  zero  on 
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the  other  edge  x  >  0,  j/  =  0.  The  resulting  wave  motion  is  a  plane  wave  at  i  =  c^t  that  carries 
a  unit  jump  in  particle  displacement  in  the  x-direction,  plus  cylindrical  dilatational  and  shear 
waves  centered  at  the  corner  of  the  quarter  plane  with  their  associated  headwaves.  A  Rayleigh 
wave  is  also  expected  on  the  free  surface  at  x  =  cut. 

For  purposes  of  the  superposition  step  indicated  in  Fig.  5  and  its  counterpart  for  the  case 
of  antisymmetric  loading,  the  only  feature  of  the  solution  of  the  quarter  plane  problem  that 
is  required  is  the  displacement  on  the  edge  y  =  0  for  all  time.  It  is  this  displacement  that 
must  be  negated  by  a  distribution  of  moving  dislocations  in  (D).  The  displacement  component 
Uy(x,0,t)  will  lead  to  the  mode  I  stress  intensity  factor,  and  likewise  the  component  Ui(x,0,  t) 
will  lead  to  the  mode  II  stress  intensity  factor. 

The  two  half  plane  problems  whose  superimposed  solutions  provide  the  solution  for  the 
quarter  plane  problem  are  labelled  (E)  and  (F);  see  Fig.  6  for  the  case  of  symmetric  loading.  In 
(E),  a  unit  step  in  normal  displacement  is  imposed  on  the  surface  at  time  t  =  0,  and  the  shear 
traction  on  this  surface  is  zero.  The  elementary  solution  is  simply  a  plane  wave  propagating  in 
the  x-direction  behind  which  u*  =  1.  This  wave  satisfies  the  quarter  plane  boundary  conditions 
on  the  edge  x  =  0,  j/  >  0  but  the  boundary  conditions  on  the  other  edge  are  violated.  The 
plane  wave  induces  no  shear  stress  cr^y  on  this  surface  but  it  does  induce  a  compressive  normal 
stress  —Xa6(t  —  x/c^)  where  S(  )  denotes  the  Dirac  delta  function.  Thus,  the  sohition  of  the 
half  plane  problem  (F)  must  satisfy  the  boundary  conditions 

<rj;y(x,0,t)  =  0,  <Tyy{x,0j)  =  XoS{t  -  xlcd) .  (2.1) 


In  considering  the  physical  dimensions  of  terms  in  (2.1)  and  of  other  fields  in  this  section,  it 
should  be  kept  in  mind  that  all  expressions  are  multiplied  by  a  unit  displacement. 

Application  of  the  Laplace  transforms  defined  in  (1.4)  to  the  wave  equations  (1.2)  and  the 
boundary  conditions  (2.1)  leads  to  the  transformed  wave  potentials  for  problem  (F)  given  by 


V{^,y,s) 


X  1  n(62  _  2^2) 
fis^  M(0 


ffiOc 


-scry 


V{^,y,s 


X  1  2n^o 

RiO 


Rc(o)  >  0 
Rc(/?)  >  0 


(2.2) 
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where  a(^)  and  i3(()  are  defined  by 


(2.3) 


and 


ff(^)  = 


1  1 
a-^'^  a  +  i 


R(0  =  (6"  -  +  4^^  a/3  . 


(2.4) 


The  superscript  F  indicates  a  feature  of  problem  (F).  The  condition  Re(a)  >  0  is  met  for  all 
points  in  the  ^-plane  by  cutting  ^-plane  along  a  <  |Re(f)I  <  oo,  Im(^)  =  0,  and  choosing  the 
branch  for  which  a  is  a  positive  real  number  at  ^  =  0.  Similar  conditions  are  imposed  on  /?. 
The  function  Jt(()  is  recognized  as  the  Rayleigh  wave  function,  and  the  only  roots  of  R(()  =  0 
in  the  complex  plane  cut  as  indicated  are  at  ^  =  ±c  =  The  domain  of  definition  of 

each  of  ^  and  ^  is  extended  from  the  common  strip  of  convergence  of  the  integral  transforms 
-a  <  Re(^)  <  a  to  the  entire  ^-plane,  cut  along  the  line  lRe(^)l  >  a,  Im(^)  =  0.  Both  ^ 
and  ijj  have  simple  poles  at  the  zeros  of  R{^).  Furthermore,  ^  has  simple  poles  at  ^  =  ±a, 
which  correspond  to  the  point  loads  on  the  surface  of  the  half  space  propagating  in  opposite 
directions  with  the  dilatational  wave  speed.  Transformed  displacements  for  (F)  can  be  obtained 
from  potentials  in  the  form 


Wx(^,0,5) 

«J^(^,0,s) 


A  1  «fr^^(2f  -  6^  +  2a[3)(M) 

M  s  R{0 

^  1  a«(0!?(C) 

/t  s  i?(0  • 


(2.5) 


It  is  noted  that  the  transformed  form  of  is  similar  to  Uy  of  Lamb’s  (1901)  problem  and  the 
transformed  form  of  is  similar  to  Uy  of  Abou-Saycd  et  al  ,  with  a  weighting  function  g{^). 

The  displacement  components  in  physical  coordinates  can  be  obtained  by  application  of 
the  Laplace  transform  inversion  integrals  and  the  Cauchy  integral  theorem.  The  displacement 
components  for  problem  (C)  are  then  obtained  by  superposition  of  the  results  from  (E)  and 
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(F),  that  is,  uf  =  ttf  +  tif .  The  result  is 


«f(x,0,0=  { 


ax  <  t  <  bx 


(2.6) 


bx  <  t  <  oo 


and 

u^(x,0,t)=-—  /  Im  I  - ^^==)  ax  <  t  <  oo  (2.7) 

^  Ja  \RiOy/a^^) 

where  integration  paths  are  taken  on  Im(^)  — ►  0“.  The  function  S{^)  is  defined  in  terms  of  the 
Rayleigh  wave  function  by 

R{0  =  K(c^-e)S{0  (2.8) 


where  k  =  2(6^  -  a?).  The  function  5(^)  has  neither  zeros  nor  poles  in  the  cut  ^-planc  and 
5(^)  -*  1  as  1^1  — ♦  oo.  With  reference  to  Fig.  5,  the  displacement  component  in  (2.6)  must 
be  negated  by  superposition  of  dislocations  in  (D).  The  method  to  be  described  in  the  next 
section  hinges  on  the  observation  that  and  are  homogeneous  functions  of  x  and  t  of 
degree  zero. 


3.  FUNDAMENTAL  MOVING  DISLOCATION  SOLUTIONS 

The  next  step  in  the  analysis  is  to  determine  the  crack  tip  stress  intensity  factor  for  a  single 
edge  dislocation  climbing  at  constant  speed  out  of  the  crack  tip;  see  problem  (D)  in  Fig.  5.  The 
body  is  stress  free  and  at  rest  up  until  the  instant  that  the  dislocation  emerges  from  the  crack 
tip.  A  solution  to  this  problem  that  is  valid  for  all  time  cannot  be  obtained  in  the  presence  of 
the  free  boundary  at  x  =  —i.  However,  a  solution  can  be  found  for  the  case  of  a  semi-infinite 
crack  in  an  otherwise  unbounded  body,  and  the  stress  intensity  factor  histories  will  be  identical 
cases  for  the  two  cases  for  time  in  the  range  0  <  t  <  2f/cd.  When  the  dislocation  begins  to 
move  from  the  crack  tip,  a  Rayleigh  wave  is  initiated  on  the  crack  face.  This  surface  wave 
travels  at  speed  cn  from  the  crack  tip  toward  the  free  edge  of  the  plate.  The  main  clfect  of  free 
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surface  in  this  case  is  to  reflect  the  Rayleigh  wave  back  toward  the  crack  tip  at  time  t  =  tjcji. 
The  reflected  Rayleigh  wave  can  have  a  significant  effect  on  the  stress  intensity  factor  when  it 
arrives  back  at  the  crack  tip.  It  is  expected  that  the  difference  in  stress  intensity  factor  histories 
for  dislocation  emission  between  the  edge  crack  problem  and  the  half  plane  crack  problem  will 
be  small  for  time  in  the  range  2//c<i  <  t  <  UfcR. 

Consider  plane  strain  deformation  of  an  unbounded  elastic  solid  containing  a  half  plane 
crack  whose  faces  are  free  of  traction.  The  material  is  initially  stress  free  and  at  rest.  At  time 
t  =  0,  an  elastic  dislocation  with  Burgers  vector  of  magnitude  2A  in  the  y-direction  begins  to 
move  at  constant  speed  v  in  the  positive  i-direction  from  the  crack  tip.  In  the  terminology 
of  elastic  dislocation  theory,  this  is  an  edge  dislocation  climbing  in  the  i-direction.  The  speed 
is  assumed  to  be  in  the  range  0  <  u  <  cj,  though  the  solution  may  be  obtained  by  the  same 
procedure  if  the  speed  exceeds  the  dilatational  wave  speed  of  the  material. 

For  the  case  of  mode  I,  the  deformation  is  symmetric  with  respect  to  the  plane  y  =  0>  Thus 
a  solution  is  sought  for  — oo  <  x  <  oo,0  <  y  <  oo,  and  0  <  /  <  oo  subject  to  the  boundary 
conditions 

ayv(x,0,<)  =  <r+(x,t) 

<riy(x,0,t)  =  0  (3.1) 


Uy{x,0,t)  =  u-{x,t)  +  AH{vt  -  x) 

where  cr+  is  the  unknown  normal  stress  distribution  on  the  plane  y  =  0  for  x  >  0,  and  =  0 
for  I  <  0.  Likewise  u_  is  the  unknown  normal  displacement  distribution  on  y  =  0  for  x  <  0, 
and  u_  =  0  for  X  >  0. 

Following  the  procedure  introduced  by  Freund  (1974),  tlie  Wiener-llopf  equation  for  this 
case  is 


E+(0  =  f/-(0  + 


A 


(3.2) 


where  E+  and  U-  are  tlie  double  Laplace  transforms  of  and  u_,  respectively,  and  d  —  \/v. 
Application  of  the  factorization  procediirc  of  the  Wiener-llojrf  method  leads  to  the  solution  for 
the  double  transform  of  <t+  in  the  form 


s+(0  =  - 


Ap«  1 


(3.3) 
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where 


and 


S*({)  =  exp|-i|^  tan-' 


62)2 


dr} 

r}±^ 


The  transform  inversion  leads  to  the  expression  for  (7+  in  the  physical  domain  as 

£r+(i,t)  =  Im  <S+(-—  4-  *0)  I  H{t  -  ax) . 

TX  1^  X  ) 


(3.4) 

(3.5) 

(3.6) 


The  stress  intensity  factor  history  for  mode  I  due  to  the  dislocation  moving  at  speed  v,  say 
ki(t;  v),  may  be  deduced  from  the  asymptotic  behavior  of  S+(^)  as  ^  ►  oo,  with  the  result 


ki(t\v) 


1 

V  r'P'  F+(d)v/t 


(A  =  l). 


(3.7) 


The  fact  that  a^{x,t)  is  a  homogeneous  function  of  its  arguments  of  order  -1  is  evident  in 
(3.6).  Furthermore,  the  dependence  of  the  stress  intensity  factor  on  the  dislocation  velocity 
V  has  been  indicated  explicitly  in  (3.7).  These  observations  are  important  in  calculating  the 
transient  stress  intensity  factor  A'/(0  for  the  crack  problem  (A).  Recall  that  this  step  requires 
the  superposition  of  dislocations  that  reproduces  a  normal  displacement  that  is  exactly  equal 
but  opposite  to  the  normal  displacement  Uy{x,0,t)  of  the  quarter-plane  problem  (C);  see  Fig.  5. 

The  case  of  mode  II  can  be  handled  in  an  analogous  way.  The  direction  of  Burgers  vector 
is  the  negative  i-direction,  rather  than  the  positive  y-direction  as  in  mode  I.  In  addition,  the 
normal  stress  (jyy  vanished  on  y  =  0  in  the  mode  II  equivalent  of  problem  (D)  instead  of  the 
component  cTxy  in  mode  I.  Suppose  that  r+(i,/)  is  tlie  unknown  sliear  stress  on  y  =  0,  .r  >  0 
in  this  case.  Without  repeating  the  details,  it  can  be  shown  that 


2'+(2^^0  =  {  77Tl\r  \f\r  r/\l  H{t-ax)  (3.8) 

rrx  b-  Us  +  «)G+(a<^-)-(fO  J  «=-«/i+i0 


where 


G+iO 


jc  +  Qs+iO 


The  fundamental  stress  intensity  factor  solution  is 


kii{t\  v) 


[2}ik  1 

V  TT  62  G+[d)sfi 


(A  =  1). 


(3.9) 


(3.10) 
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The  fundamental  stress  intensity  factor  solutions  (3.7)  and  (3.10)  vvill  be  used  in  the  next  section 
to  construct  the  full  mode  I  and  mode  II  stress  intensity  factor  histories. 

4.  MIXED  MODE  STRESS  INTENSITY  FACTORS 

In  the  preceding  section,  the  stress  intensity  factor  due  to  the  motion  of  a  dislocation  away 
from  the  tip  of  the  crack  at  constant  speed  was  considered.  Suppose  now  that  a  continuous 
distribution  of  dislocations  is  emitted  from  the  crack  tip,  each  moving  with  a  different  constant 
speed.  Furthermore,  suppose  that  the  time  of  emission  and  the  amplitude  of  the  Burgers 
displacement  both  depend  on  the  speed  of  the  dislocation.  For  any  continuous  variation  of 
the  amplitude,  say  ‘w(v),  and  the  time  of  emission,  say  to(v),  with  speed  v  the  continuous 
distribution  represents  a  displacement  distribution  that  is  a  homogeneous  function  of  degree 
zero  of  position  and  time.  From  (3.7)  or  (3.10),  the  stress  intensity  factor  resulting  from  this 
dislocation  distribution  is 

K(t)=  -  1^'  k{l-t,-,v)^^dv  (4.1) 

where  ui  and  V2  are  appropriate  limits  on  the  actual  range  of  speed. 

In  the  present  case,  the  distribution  of  climbing  edge  dislocations  giving  rise  to  an  opening 
displacement  across  the  plane  y  =  0  results  in  the  mode  I  stress  intensity  factor.  Thus,  to  in 
(4.1)  is  replaced  by  Uy  from  (2.7).  If  the  boundary  loading  acts  at  time  t  —  0  and  at  place 
X  =  then  the  time  at  which  the  displacement  level  with  speed  x/t  arrives  at  i  =  0  is 
to{v)  =  (1  —  u/c<f)f/u.  The  first  displacement  level  to  arrive  at  x  =  0  does  so  with  speed  Cd, 
and  the  range  of  speeds  of  all  displacement  levels  involved  at  time  t  is  C/t  <  r  <  Cd-  Therefore, 
the  mode  I  stress  intensity  factor  resulting  from  negation  of  the  displacement  distribution  , 
and  the  corresponding  mode  II  stress  intensity  factor  resulting  from  negation  of  ,  are 

(4.2) 

K,i{t)  =  ~£  (4.3) 

where  h"  —  t/t  and  the  inverse  speed  h  =  \fv  has  been  used  as  the  integration  variable 
for  convenience.  The  integral  expression  for  A'/  can  be  evaluated  by  using  the  fundamental 
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formula  for  A:/(^;  v)  in  (3.7)  and  the  explicit  expression  for  the  displacement  Uy  given  in  (2.7). 
The  procedure  is  described  by  Freund  (1974),  and  only  the  result  is  included  here,  that  is. 


Ki{t)  = 


K  f h*  2a^(2h^~b^)P{h)  .. 

(h+a)^{h-a)(lr-h) 


— 
yJ'l'Kl  I 


^/7a 


+ 


^/i*+aS+(a)  c+o  Ti^S— (c)(c^  — a*) 


_ 


a  <  /i*  <  6 


b  <  h*  <  c 
c  <  h*  <3a 


(4.4) 


The  history  of  Kj  is  square  root  singular  at  t/£  =  c,  but  a  nonzero,  finite  value  is  obtained  at 
t/(.  =  a  from  the  integration.  Following  a  similar  process,  Ku  is  obtained  as 


/  K  rh-  4^hP{h)./b^ 

I  *■  v'(Ai2-a2)(A*-A) 


a  <  h*  <  b 


Knit)  =  - 


A 

V2iri 


—  f  ^  dfl 

^  y/(h^-a^)(h‘-h) 

—  f** 

^  •*“  y(h2-a2)(/‘*-'‘) 


b  <  h’  <  c 

■  (4.5) 


I  2o^c  2c^-6^-2v/^(c»-a»)(e^-6^)  i 

■^■6^(02-02)  P+(c)S-(c)  x/h--c 


c  <  h*  <  3a 


The  history  of  Ku  has  a  square  root  singularity  at  tji  -  c,  but  again  a  finite  value  is  found  at 
tjt  =  a  from  the  integration.  In  the  above  e.xpressions,  F(/i)  =  (c  +  li)(2lr  -  b-)S+(h)/[i2h?  - 

62)4  +  16/l<(h2  _a2)(i2 

For  numerical  evaluation  of  these  c,xpressions  Poisson’s  ratio  i/  is  assumed  to  be  equal  to 
0.25.  In  this  case,  the  ratios  of  the  slowness  arc  b~  =  3a^  and  =  3.519  a'.  The  stress  intensity 
factors  have  been  evaluated  numerically  by  application  of  the  appropriate  Ga.ussian  quadrature 
rules  with  Chebychev  polynomial  interpolations.  Tiie  results  are  shown  in  Fig.  7  and  Fig.  8. 
The  results  are  normalized  by  a  factor  that  is  the  static  mode  11  stress  intensity  factor  as  t  -+  oo 
for  the  problem  posed  in  section  1  and  shown  schematically  in  Fig.  3.  This  stress  intensity  factor 
is  easily  calculated  to  be  F/{2(1  -  p^)\/TrI]  (times  the  unit  imposed  displacement)  by  means  of 
the  Af-integral  procedure  introduced  by  Freund  (1978). 

Up  to  this  point,  the  stre.ss  pul.se  grazing  the  lower  edge  of  the  crack  surface  has  had  the 
form  of  a  displacement  with  step  function  time  dependence.  The  stress  intensity  factors  for 
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general  time  dependence  can  be  obtained  by  convolution.  If  tlic  velocity  on  the  lower  edge  is 
prescribed  as  an  arbitrary  function  of  time  then  the  stress  intensity  factor  for  this  situation 
is 

AC(t)=  /  K{t-s)v{s)ds.  (4.6) 

JQ 

For  special  example,  suppose  that  the  velocity  is  a  step  function  v{t)  =  Von{t),  which  is 
essentially  the  loading  condition  used  by  KalthofF  (1987).  The  resultant  values  of  /C/(t)  and 
^//(O  obtained  according  to  (4.6)  are  presented  in  Fig.  9  and  Fig.  10.  The  normalizing  factor 
here  is  Evoy/t/ir  /[2cd(l  -  i/*)]. 

5.  DISCUSSION 

In  the  foregoing  sections,  the  exact  elastodynamic  stress  intensity  factor  history  has  been 
determined  for  a  particular  case  of  asymmetric  impact  of  an  edge  cracked  plate.  The  stress 
intensity  factor  results  for  impulsive  imposed  velocity  v{t)  =  6{t)  are  presented  in  Figs.  7  and 
8.  In  themselves,  these  results  are  of  limited  practical  significance.  However,  they  provide 
the  fundamental  building  block  for  determining  the  transient  stress  intensity  factors  for  more 
realistic  imposed  ve’ocity  boundary  conditions.  In  addition,  the  existence  of  an  exact  solution 
for  a  configuration  of  this  sort  provides  a  valuable  check  on  numerical  procedures  that  are  being 
developed  to  analyze  more  complex  geometries  and,  eventually,  impact  with  nonlinear  material 
response. 

The  exact  stress  intensity  factor  histories  for  step  velocity  loading  v(t)  =  Vo  II {t)  are  given 
in  Figs.  9  and  10.  As  anticipated,  when  the  imposed  velocity  on  the  boundary  is  in  a  direction 
into  the  material,  then  the  adjacent  crack  face  bulges  outwaid  and  the  resulting  mode  I  stress 
intensity  factor  is  negative.  If  the  crack  is  initially  closed  the  faces  would  press  against  each 
other  and,  in  fact,  no  mode  I  stress  intensity  factor  would  develop.  On  the  other  hand,  if  the 
crack  is  initially  held  open  by  some  background  equilibrium  load  or  if  the  geometry  is  such  that 
a  small  opening  exists  then  the  crack  will  indeed  tend  to  close  upon  application  of  the  boundary 
velocity. 

The  magnitude  of  the  mode  I  stress  intensity  factor  is  always  less  than  the  magnitude  of  the 
corresponding  mode  II  stress  intensity  factor,  but  the  value  is  significant.  A  single  parameter 
that  is  useful  in  characterizing  the  near  tip  field  under  mixed  mode  conditions  was  introduced 
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by  Shih  (1974).  The  parameter  is  defined  by 


and  it  is  called  the  mixity  parameter.  The  dependence  of  the  mixity  parameter  on  time  for  the 
problem  analyzed  here  is  shown  in  Fig.  11  where  it  is  seen  that  the  value  varies  only  a  small 
amount  from  0.25  over  the  time  interval  of  interest. 

If  it  is  assumed  that  the  crack  will  grow  in  a  direction  determined  by  the  condition  that  the 
circumferential  tensile  stress  within  the  asymptotic  field  is  maximum  then  the  angle  between 
the  crack  line  and  the  direction  of  growth  satisfies 

3sin^(^)  -  tan(|A4')cos(^)sin(^)  -1  =  0.  (5.2) 

The  angle  0e  that  satisfies  this  condition  for  the  problem  analyzed  here  is  also  shown  in  Fig. 
11  where  it  is  seen  that  this  angle  varies  only  slightly  from  63°  in  the  time  interval  of  interest. 

The  possible  influence  of  crack  face  interaction  resulting  from  the  tendency  for  the  crack  to 
close  under  compressive  loading  is  difficult  to  assess  at  this  point.  In  the  pure  mode  I  problem 
considered  by  Freund  (1974)  the  crack  could  be  view  as  closed  with  normal  pressure  acting  or 
open  with  no  traction  acting.  In  the  mixed  mode  case  considered  here,  however,  the  crack  faces 
would  tend  to  slide  with  respect  to  each  other  if  they  are  in  contact  so  that  the  two  modes  are 
coupled  in  a  way  that  has  not  yet  been  sorted  out.  Perhaps  some  progress  on  this  point  could 
be  made  by  considering  the  case  of  closed  crack  faces,  that  is,  continuous  normal  traction  and 
velocity  across  the  crack,  but  zero  shear  traction  over  the  entire  crack  surface.  This  problem 
could  be  analyzed  by  the  procedure  used  here. 

Data  taken  from  the  work  of  Kalthoff  (1987)  is  also  shown  in  Fig.  10.  He  measured  the 
speed  of  the  impacting  projectile  in  the  configuration  of  Fig.  1,  and  the  resulting  stress  intensity 
factor  histories  wore  observed  by  means  of  the  optical  shadow  spot  method.  In  his  experiments, 
interaction  between  the  crack  faces  was  avoided  by  cutting  the  cracks  into  the  plate  so  that 
there  was  a  finite  gap  between  the  crack  faces.  Thus,  both  a  mode  II  stress  intensity  factor 
and  a  compressive  mode  I  stress  intensity  factor  could  be  generated,  as  modeled  in  the  present 
analysis.  To  enter  this  data  for  steel  specimens  struck  by  steel  projectiles  on  Fig.  10,  the  velocity 
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value  Vo  imposed  on  the  edge  of  the  specimen  was  determined  by  assuming  that  the  portion 
of  the  specimen  struck  by  the  projectile  (see  Fig.  1)  had  the  same  elastic  impedance  as  the 
projectile  itself.  Thus,  from  elementary  wave  propagation  theory,  the  value  of  Vo  was  essentially 
one-half  of  the  speed  of  the  incident  projectile.  The  agreement  between  the  analytical  model 
results  and  the  experimental  results  for  the  low  velocity  result  {vo  =  6.5  m/ s)  is  excellent  over 
the  time  range  for  which  data  are  reported.  On  the  other  hand,  the  observed  mode  II  stress 
intensity  factor  for  the  higher  velocity  impact  (vo  =  16.5  m/s)  is  less  than  the  theoretical  result 
up  until  the  time  Cdt/t  ~  2,  and  substantially  above  thereafter.  Kalthoff  (1987)  also  reported 
results  of  one  experiment  on  the  brittle  plastic  Araldite  B  which,  when  normalized,  were  closer 
to  the  analytical  results  than  the  high  impact  velocity  steel  results  but  which  differed  in  the 
same  general  way.  Nonetheless  the  trend  between  the  observed  variation  of  Kn  with  time  and 
the  analytical  result  are  consistent  in  both  cases  up  until  Cit/l  ~  2.  The  differences  thereafter 
could  be  due  to  nonlinear  material  response  in  the  crack  tip  region  or  to  the  influence  of  other 
features  of  the  experiment  not  taken  into  account  in  the  analytical  model,  such  as  the  presence 
of  the  second  crack.  Both  of  these  possibilities  arc  being  pursued  through  further  work  based 
on  detailed  numerical  simulation  of  the  process. 

Some  understanding  of  the  role  of  crack  tip  processes  in  the  interpretation  of  dynamic 
fracture  data  obtained  by  means  of  the  optical  shadow  spot  method  has  been  achieved.  The 
shadow  spot  image  on  a  screen  or  camera  is  created  as  a  locus  of  points  of  multiple  mapping 
of  a  family  of  light  rays  influence  by  deformation  of  the  specimen.  The  curve  on  the  specimen 
that  is  geometrically  mapped  into  the  boundary  of  the  shadow  spot  must  be  within  a  region 
where  the  mechanical  fields  are  adequately  represented  by  the  asymptotic  elastic  crack  tip  field, 
and  this  condition  imposes  quite  severe  constraints  on  the  method.  The  presence  of  a  crack 
tip  plastic  zone  can  distort  the  shadow  spot  and  increase  its  size  so  that  the  inferred  value  of 
stress  intensity  factor  is  not  accurate  (Rosakis  and  Freund,  1981).  In  the  experiments  reported 
by  Kalthoff  (1987)  a  well-defined  shear  band  grew  from  the  crack  tip  in  tlie  plane  of  the  crack 
for  the  case  of  high  velocity  impact,  and  it  is  possible  that  the  presence  of  the  band  influenced 
the  inferred  stress  intensity  factor  values. 
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FIGURE  CAPTIONS 


Figure  1. 

Figure  2. 

Figure  3. 

Figure  4. 

Figure  5. 

Figure  6. 

Figure  7. 

Figure  8. 
Figure  9. 

Figure  10. 

Figure  11. 


Experimental  configuration  introduced  by  Kalthoff  and  Winkler  (1987)  for  dy¬ 
namic  fracture  initiation  studies  under  mixed  mode  conditions. 

Pattern  of  wave  fronts  generated  in  the  portion  of  the  specimen  shown  in  Fig. 
1  between  the  two  cracks  due  to  normal  impact  of  the  projectile. 

A  schematic  representation  of  the  boundary  value  problem  considered  here, 
showing  the  edge  crack  of  length  I  and  the  imposed  normal  velocity  imposed 
asymmetrically  on  the  edge  of  the  plate. 

Decomposition  of  the  asymmetric  boundary  value  problem  in  Fig.  1  into  a 
symmetric  loading  system  (A)  and  an  anti-symmetric  loading  system  (B). 
Decomposition  of  problem  (A)  into  a  quarter  plane  problem  (C)  and  a  dynamic 
dislocation  propagation  problem  (D). 

Decomposition  of  the  quarter  plane  problem  (C)  into  two  half  plane  problems, 
each  of  which  can  be  solved  by  standard  methods  of  analysis. 

The  mode  I  stress  intensity  factor  for  the  case  of  impulsive  imposed  velocity 
v(<)  =  6{t),  normalized  by  E‘/[2(l  -  versus  normalized  time. 

The  mode  II  stress  intensity  factor  corresponding  to  the  result  in  Fig.  7. 

The  mode  I  stress  intensity  factor  for  the  case  of  step  loading  v{t)  =  Vq  H{t), 
normalized  by  Evoy/Tj^ /[2(1  -  versus  normalized  time. 

The  mode  II  stress  intensity  factor  corresponding  to  the  result  in  Fig.  9.  The 
data  points  are  those  reported  by  Kalthoff  (1987). 

The  mixity  parameter  defined  in  (5.1)  for  the  stress  intensity  factors  given  in 
Figs.  9  and  10,  and  the  corresponding  direction  of  crack  growth  based  on  a 
maximum  circumferential  tensile  stress  criterion. 
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